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Abstract 

We derive Gamow vectors from S-matrix poles of higher multiplicity in analogy to 
the Gamow vectors describing resonances from first-order poles. With these vectors 
we construct a density operator that describes resonances associated with higher order 
poles that obey an exponential decay law. It turns out that this operator formed by 
these higher order Gamow vectors has a unique structure. 
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1 Introduction 



Resonances in quantum mechanics can be described by poles of the analytically continued 
S-matrix on the second sheet of a two-sheeted Riemann surface (Newton, 1982; Goldberger 
and Watson, 1964; Bohm 1979, 1980, 1981, 1993). These poles appear in conjugate pairs 
below (at zr = Er — iT/2 ) and above (at z* R = Er + iT/2) the positive real axis, where the 
pole at zr corresponds to a decaying state at times t > and the pole at z* R to the growing 
state at t < 0. These poles lead to the description of the Gamow vectors with energy 
E R and lifetime r = h/T. The Gamow vectors possess all the properties of resonances, 
in particular, an exponential decay law and a Breit-Wigner energy distribution. One can 
extend this derivation of the Gamow vectors from first order poles to S-matrix poles of higher 
multiplicity (Antoniou and Gadella, 1995; Bohm et al, 1997). An S-matrix pole of order r 
at the position zr = Er — iT/2 on the second Riemann sheet leads to a set of r generalized 
eigenvectors of the Hamiltonian of order k — 0, 1,..., r — 1, which are Jordan vectors of 
degree k + 1 to the generalized eigenvalue E R — iT/2 and which are elements of a generalized 
complex eigenvector expansion (nuclear spectral theorem in the rigged Hilbert space). The 
form of this generalized complex eigenvector expansion suggests the definition of a state 
operator (density matrix) for the microphysical decaying state from a higher order pole. This 
microphysical state is not a pure state, but a mixture of non-reducible components. In spite 
of the fact that the k-th. order Gamow- Jordan vectors have a polynomial time-dependence 
besides the exponential, which in the past were always associated with resonances from 
higher order poles, this microphysical state obeys a purely exponential decay law. 

Resonances from higher order poles, in particular double poles, have already been de- 
scribed about 30 years ago (Goldberger and Watson; Newton, 1982; Goldhaber, 1968), but 
were always associated with an additional polynomial time dependence not confirmed in 
experiment. However, operators containing finite dimensional matrices consisting of non- 
diagonalizable Jordan blocks have been discussed in connection with resonances numerous 
times in the past (Mondragon, 1994; Stodolsky, 1970; Lukierski, 1967; Dothan and Horn, 
1970; Katznelson, 1980; Bhamathi and Sudarshan, 1995; Brandas and Dreismann, 1987; 
Antoniou and Tasaki, 1993). Using the formalism of the rigged Hilbert space, I. Antoniou 
and M. Gadella derived the Gamow- Jordan vectors, or higher order Gamow vectors, from 
the higher order poles of the S-matrix. 

In the second section of this article, we recall some of the notation (Bohm, 1993; Bohm et 
al, 1997) needed for the description of the scattering experiment in the rigged Hilbert space 
formalism. We show how to obtain these hypothetical vectors associated with the higher 
order S-matrix poles and show that they are Jordan vectors (Baumgartel, 1984; Kato, 1966; 
Gantmacher, 1959). We derive their properties under the action of the Hamiltonian and their 
semigroup time evolution. In the third section, we discuss possible operators formed by these 
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vectors to represent microphysical states describing resonances. In the fourth section, we 
will ask for the converse : Going out from an exponential decay law for the time evolution of 
resonances from higher order S-matrix poles, what is the most general form of the operator 
formed by dyadic products of higher order Gamow vectors? 

2 Poles of the S-Matrix and Gamow-Jordan Vectors 

We recall that in the rigged Hilbert space formalism one uses two different space triplets for 
the set of in-states <p + G C 7i C describing the preparation process of the scattering 
experiment and the set of observables G $ + C Ji C describing the registration 
process (Bohm et al, 1997; Bohm and Gadella, 1989) The in-state + , which evolves from 
the prepared in-state m outside the interaction region, is determined by the accelerator. 
The so-called out-state ip~ (or ^ ont ) is determined by the detector. \^ ont ) (^ out \ is therefore 
the observable which the detector registers and not a state. The S-matrix elements are given 
by the projection of the set of out-states {0 out } onto the set of observables {i/j ont } (Bohm, 
1993) 

(^ out ,0 out ) = (^ ou \S0 in ) = = / dE (^\E-)S(E)( + E\(f) + ) . (2.1) 

J spec H 

The vectors |i? ± ) G $ x are the scattering states (Dirac kets) and are eigenvectors of the 
exact Hamiltonian with energy label E, which can take values on a two-sheeted Riemann 
surface. We choose to ignore all other labels of the basis vectors {E^), since nothing im- 
portant is gained in our discussion if we retain the additional quantum numbers, e.g., the 
angular momentum quantum numbers / and I3 or the polarization or channel quantum num- 
bers 77. Thus, we shall restrict our discussion on one initial channel r] — tja and one final 
channel 7/ = t\b (e.g., t)b = t]a for elastic scattering), and we shall consider the Z-th partial 
wave of the 7773-th channel (Bohm, 1993): i.e., S(E) = S T [ ,B (E). We consider the model in 
which the S-matrix is analytically continued to a two-sheeted Riemann surface in the energy 
representation (Bohm, 1993) and in which the S-matrix S(u), u G(E, has one r-th order pole 
at the position uo = z R (z R = E R — iT/2) in the lower half-plane of the second sheet, (and 
consequently there is also one r-th order pole in the upper half-plane of the second sheet at 
uo = z* R ). In this paper we will not discuss the pole at z R , as it leads to r growing higher 
order Gamow vectors, and the correspondence between the growing and decaying vectors is 
just the same as for first-order pole resonances (r = 1). The model that we discuss here can 
easily be extended to any finite number of finite order poles in the second sheet below the 
positive real axis. 
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The unitary S-matrix of a quasistationary state associated with an r-th order pole at 
z R = E R — iT/2 in the lower half-plane of the second sheet (denoted by II) is represented by 
(Bohm, 1993, sect. XVIII.6) 

Su(cu) = e 2 ^) e 2i7(aj) , (2.2) 

where S R (cu) = r arctan( 2( - £ ^_^ ) is the rapidly varying resonant part of the phase shift, 
and 7(0?) is the background phase shift, which is a slowly varying function of the complex 
energy to. r is a dimensionless quantity that, due to the analyticity properties of the S- 
matrix (Bohm, 1993, section XVIII.6), takes integer values, where r > leads to a decaying 
resonance of order r, and r < to its corresponding growing state. Using the identity 
arctan-g^-^j = | (ln(co> — E R — iT/2) — \n(uj — E R + iT/2)), one rewrite Sn(u): 

Su (lu) - ( — - Er - lV/2 \ T c 2 ^ ] - e 2i ^ c 2 ^ (2 1) 

" \u-(E R -iT/2)) 6 ~ e {u-z R y e (2 - 3) 



We insert this into ( |2.1|) and deform the contour of integration C_ through the cut along the 
spectrum of H into the second sheet (Bohm, 1979, 1980, 1981, 1993). Then one obtains 

(<T,0 + ) = [ duo {ij-\uo-)S ll {uo)( + uo\ ( t> + )+ (2.4) 



+ E ( n + 1) (~ iv r +1 £ {u _ ZR)n+1 <w>; m«o < o . 

The first integral does not depend on the pole and is called the "background term". The 
contour C_ can be deformed into the negative axis of the second sheet from to — oon. 



— OO0 

k+ 



(iP~, <P + ) = / dE(4,-\E-) Sn(E) ( + E\<P + ) + (tP~, + )p. t . (2.5) 



higher order pole term (ip , 4> + )p.t. we obtain using the Cauchy integral formulas §^ rjz^m+i duJ 
2 ^r f {n \z)\ z=ZR where /<»)(*) = 

(^-,0 + )p.t. = E(„ + i) (~^ n+1 (~ 2 -f) e 2l7H (W))^ (2-6) 

n=0 ^ ' \ ■ J u ZR 



where (. . . )^/l ZR denotes the n-th derivative with respect to lu taken at the value u = z R . 



Since the kets \uj~) are (like the Dirac kets \E~)) only defined up to an arbitrary factor or, 
if their "normalization" is already fixed, up to a phase factor, we can absorb the background 
phase e 2 * 7 ^ into the kets \uj~) and define 

\u f ) = |cj-)e 2i7(w) . (2.7) 

Note that this phase is not trivial, e.g. \E+) = \E~)Sn(E) = \E~)e 2iSR ^e 2i ^ E \ except for 
the case when the slowly varying background phase 7(0;) is constant and the \uj 1 ) are identical 
with \uj~) up to a totally trivial constant phase factor. In general (|2.7|) is a non-trivial gauge 
transformation. We will keep the phase 7 in our Gamow vectors, but not investigate their 
properties further. The results of this section are true also for the case when |4) are 
exchanged by \z~j^ k \ But as we cannot just ignore the existence of the background integral 
in (2.4), one has to keep in mind that their existence is not irrelevant, if one deals with poles 
of order r > 1. 

Taking the derivatives, we write the pole term as 

= E(„+i) h t >" +1 (-f?) E © <*-i*#* ,( "^ + *«i* + >. (2-8) 

Tl — k — 

where we denote the n-th derivative of the analytic function {il)~\z~ 1 ) by {il)~\z 1 ) (jl " > with 
value ( , - l4) at z = zr. Since (ip~~\E~) E S R 7Y 2 , i.e. element of the Schwartz space 
and of Hardy class (Duren, 1970; Hoffman, 1962; Bohm and Gadella, 1989), it follows that 
{il)~\z 1 )^ is also an analytic function in the lower half-plane of the second sheet, whose 
boundary value on the positive real axis (ip~ ^E 1 )^ E S fl 7Y 2 . Analogously, we denote by 
( n \ + z\cf) + ) the n-th derivative of the analytic function ( + z\(f) + ). Again, ^ n \ + z\(f) + ) is analytic 
in the lower half-plane with its boundary value on the real axis being ( n \ + E\(f) + ) E S fl 7i 2 _. 
The r-th order pole is therefore associated with the set of r generalized vectors 

I4) (0) , I4> (1) , •••,l4f ) ,---,l4) (n) - (2.9) 

For the first order resonance pole this, of course, reduces to the single vector \z^} = I4) < "°' ) 
in agreement with (Bohm, 1979, 1980, 1981, 1993). 

We can now establish the complex basis vector expansion in analogy to the Dirac ba- 
sis vector expansion (Nuclear Spectral Theorem in the rigged Hilbert space, (Gel'fand and 
Vilenkin, 1964; Bohm and Gadella, 1989)). If we return to the complete S-matrix ele- 
ment ( |2.4j ) and insert the pole term, we get 

/•-oon 

(ip-,(f> + ) = / dE(*l)-\E + ){ + E\(t) + ) (2.10) 







-EC+i) H-(- ir >"E it) ( "-"<w> 

n=0 ^ ' ' k=0 ^ ' 
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Omitting the arbitrary if) G $+ and rearranging the sums in the second term, we obtain 
the complex basis vector expansion for an arbitrary <p + G 

0+ = / dE\E + ){ + E\cf) + ) + £ , (2.11) 



where the coefficients are given by 

r-1 



n=k 



r \ n 
n+l) [k 



k=0 



L{n - k \ + z R \^ + ) (2.12) 



This complex generalized basis vector expansion is the most important result of our ir- 
reversible quantum theory (as is the Dirac basis vector expansion for reversible quantum 
mechanics). It shows that the generalized vectors (|2.9| ) (functionals over the space $ + ) are 
part of a basis system for the + G $_ and form together with the kets \E + ), — oon < E < 0, 
a complete basis system. The vectors (|2.9| ) span a linear subspace Ai ZR C of dimension 
r: 



£ = £ai4> (fc) ; ae(D (2.13) 



A;=0 



If there are iV poles at zr { of order rj, then for every pole there exists a linear subspace 

Note that the label k of the higher order Gamow vectors is not a quantum number in 
the usual sense. Basis vectors are usually labeled by quantum numbers associated with 
eigenvalues of a complete system of commuting observables (Bohm, 1993, chap. IV). But 
there is no physical observable that the label k is connected to. Therefore, the different 
14) in the subspace M. ZR do not have a separate physical meaning. 

Now that (|2.11|) has established the generalized vectors ( |2.9|) as members of a basis 
system (together with the \E + ); — oon < E < 0) in <3>^, we can obtain the action of the 
Hamiltonian H by the action of the operator H x on these basis vectors. We will write this 
Hamiltonian in terms of its matrix elements in this basis. For this purpose we replace the 
arbitrary ip~ G $ + in ( 2.10Q by if)~ = Hip~ which is again an element of $ + , and one finds 



(Antoniou and Gadella, 1995; Bohm et al, 1997) 

= ^(^-|4) (fc) + ^14) (fe_1) ; fc = i,...,r- 

(<ri4> (0) = wi4> (0) = (if)-\H x \zip = zniriziP ■ 



(2.14) 
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This can also be written as a set of functional equations over $_ 



'-R 

H X W R ) 



z R \zlf) + k\zlf-V 
z rWr) ■ 



k 



(2.15) 



This means that H x restricted to the subspace M. ZR is a Jordan operator of degree r, and the 
vectors |4) > k = 0, 1, 2, . . . , r — 1 are Jordan vectors of degree k + 1 (Baumgartel, 1984; 
Kato, 1966; Gantmacher, 1959). They fulfill the generalized eigenvector equation (Lancaster 
and Tismenetsky, 1985), 



(ij x -^) fe+1 l4f) = o. 



Since, according to (|2.11|) , the basis system also includes the \E + ), 
indicate this by a continuously infinite diagonal matrix equation 



(Vvi£ + >) = ((r\H\E+)) = (V)(Vi£ + >) 



(2.16) 

oo n < E < 0, we 
(2.17) 



where (^(ip \E + )^j indicates a continuously infinite column matrix. Then ( [2.14j ) can be 



rewritten as: 



h- 14 



)(o) x 



z D p 



(^j-\H x \z"< R f 
\ (^~\H X \E+) ) 



fz R 

1 Z R 

2 


V0 






Zr 








o \ / y>- 



1 z R 
(£)/ 



I4f) \ 



)(D 



(^-I4) (r - 1} 

V <ri£+> / 



(2.18) 



In this matrix representation of H x , the upper left r x r submatrix associated with the 
complex eigenvalue zr is a (lower) Jordan block of degree r. One can attain the standard 
form of a Jordan block with l's on the lower diagonal by simply choosing the normalization 



I4) (fc) - M l4) (fc) 



and 



(0/h 



(0/- 



. 1 



(2.19) 



Next, we are going to discuss the time evolution of the higher order Gamow vectors. We 



replace the arbitrary ip G $ + in Q2.10 ) by ip = e lHt ip . We recall that e lHt needs to be 



a continuous operator with respect to the topology r$ + of the space $ + , and its values e lzt 
need to be holomorphic in all $ + . Its conjugate (e lHt ) x which acts on the vectors |u;~} G $ + 
is only defined for positive values of the parameter t (semigroup time evolution), 



(e iHt ip~\cj~) = (e im ) 



-iujt 



t > 



(2.20) 
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for all \oj~) G $*. Then 

k 

iz R t 



e 

p=0 



E (^) HW14) M ■ (2-21) 

p=0 

This can be written functional equation as 

(e***) X |4) W = ^ E (!) H*)*"* I4) (P) • (2-22) 



p 

In the same way one derives for the complex conjugate: 

Wpznle™ = e iz ^ E (') (^~ 9 (q) Cz R \ . (2.23) 

The same formulas apply also to the vectors \z^f k ^ and ( 1 \~zr\ with background phase 7 = 0. 
It is important to note that the time evolution operator {e tHt } x transforms between different 
14) j or different \z~j^f k \ k = 0, 1 . . . , n, that belong to the same pole of order r at z = zr, 
but the time evolution does not transform out of M. ZR - 



3 States from Higher Order Gamow Vectors 

Gamow states of zeroth order with their empirically well-established properties (exponential 
time evolution, Breit-Wigner energy distribution) have been abundantly observed in nature 
as resonances and decaying states. Theoretically, there should be no reason why quasista- 
tionary states (i.e. states that also cause large time delay in a scattering process (Bohm, 
1993, chap. 18) associated with integers r > 1 in ( |2.3| ) should not exist. However, no such 
quasistationary states have so far been established empirically. One argument against their 
existence was that the polynomial time dependence, which was always vaguely associated 
with higher order poles (Goldberger and Watson; Newton, 1982; Goldhaber, 1968), and 
which has not been observed for quasistationary states. The question that we want to dis- 
cuss in this section is, whether there is an analogous physical interpretation for the higher 
order Gamow vectors as for the ordinary Gamow vectors, namely as states which decay (for 
t > 0) or grow (for t < 0) in one prefered direction of time ("arrow of time") and obey the 
exponential law. 
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In analogy to von Neumann's definition of a pure stationary state using dyadic products 
\f){f\ of the energy eigenvectors |/) in Hilbert space, microphysical Gamow states connected 
with first order poles can be defined as dyadic products of zeroth order Gamow vectors 
(Bohm, 1979, 1980, 1981, 1993; Bohm et al, 1997), 

W^ = \z R )(-z R \. (3.1) 

(Note that this definition is not related to the scattering background phase 7 which entered 
in (|2.7D . ) The time evolution of this Gamow state is exponential, 

W G (t) = {e im ) x \z- R ){-z R \e im 

= e-^lz^i'znle^ (3.2) 

= e~ i{ER ~ i{r/2))t \z R }(~z R \ e^fl+^r/ 2 ))* 

= e- rt W G (0) } t>0. 

Mathematically, equation (|3.2j ) is understood as the functional equation of 

(ip-\W G (t)\ilj-} = e~ rt {ijj-\W G \^-) for all ip~ G $+ and t > 0. 

This shows how important it is in the RHS formalism of quantum mechanics to know what 
question one wants to ask about a microphysical state when one makes the hypothesis (|3.1|) . 
The vectors ip~ G represent observables defined by the detector (registration apparatus), 
and therefore the operator W G represents the microsystem that affects the detector. There- 
fore the quantity (ip~\W G \ip~) is the answer to the question: What is the probability that the 
microsystem affects the detector? 

If the detector is triggered at a later time t, i.e. when the observable has been time translated 

hrXtM — e im \r){r\e- mt = \r(t))(r(t)\ (3.3) 

then the same question for t > has the answer: The probability that the microsystem 
affects the detector at t > is 

(*lj-(t)\W G \tlj-(t)} = (e iHt ip-\W G \e iHt ip-) 

= {i>~\ (e iHt ) x W G e iHt \^-) (3.4) 
= e- rt (ip-\W G \tlj-} . 

This means that (|3.4j) is the probability to observe the decaying microstate at a time t 
relative to the probability (ip~\W G \ip~) at t — (which one can "normalize" to unity by 
choosing the appropriate factor on the right-hand side of ( |3.1|) ). 
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The question that one asks in the scattering experiment is different. There the pole term 
(P.T.) of (|2.5| ) for r = 1 describes how the microsystem propagates the effect which the 
preparation apparatus (accelerator, described by the state <ft + ) causes on the registration 
apparatus (detector, described by the observable This involves both the observables 

ip~ G $+ and the prepared states 4> + G and one would ask the question: What is the 
probability to observe ip~(t) in a microphysical resonance state of a scattering experiment 
with the prepared in-state (p + ? 

In distinction to the decay experiment where one just asks for the probability of ip~ G 
in the resonance scattering experiment one asks for the probability that relates ip~ G $+ to 
+ G $_ via the microphysical resonance state. Therefore the mathematical quantity that 
describes the microphysical resonance state in a scattering experiment cannot be given by 
\ z r)(~ z r\i but must be given by something like |zr)( + 2r|- The probability to observe ip~ in 
the prepared state <j) + , independently of how the effect of <fr + is carried to the detector ip~ , is 
given by the S-matrix element ( |2.1| ), 4> + )\ 2 - The probability amplitude that this effect 
is carried by the microphysical resonance state is then given by its pole term , + )p.t.- 
In analogy to the decay experiment one can also compare these probabilities at different 
times. For this purpose one translates the observable ip~ in time by an amount t > 0, 

V>~ — > ip~(t) = e iHt ip~; t>0 (3.5) 

Physically, this would correspond to turning on the detector at a time t > later than for 
ip~ . One obtains 

(r(*),0 + ) P . T . = -27rr(e^Vl4X + ^i^ + ) 

= -27rre-^<^l4>( + ^|0 + > (3.6) 
= e-^e- r */ 2 (^,0 + ) P . T . 

This means that the time dependent probability, due to the first order pole term, to measure 
the observable tp~{t) in the state <p + is given by the exponential law: 

|(e^V-,0 + )p.T.r = e- r '|(r,0 + )p.T.r. (3.7) 

This is as one would expect it if the action of the preparation apparatus on the registration 
apparatus is carried by an exponentially decaying microsystem (resonance) described by a 
Gamow vector. 

Therewith we have seen that there are two ways in which a resonance associated with a 
first order pole of the S-matrix (r = 1) can appear in experiments, and therefore there are 
two different forms to represent the decaying Gamow state: 

in a decay experiment : \ z r){ ~ z r\ (3.8a) 

and in a scattering experiment : (3.8b) 
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An analogous statement holds for the Gamow states associated with the pole in the upper 
half-plane. The first representation is the one used in the S-matrix when one calculates the 
cross section; the second representation is the one used when one calculates the Golden Rule 
(decay rate). In contrast to von Neumann's formulation where a given state (representing an 
ensemble prepared by the preparation apparatus) is always described by one and the same 
density operator |/)(/|, the representation of the microphysical state in the RHS formulation 
depends upon the kind of experiment one performs. That a theory of the microsystems must 
include the methods of the experiments has previously been emphasized by G. Ludwig. 

We will now discuss the mathematical representations of Gamow states associated with 
higher order poles of the S-matrix (r > 1) under the two aspects above. 
In analogy to the case for r = 1, we conjecture that the representation for the microphysical 
system in the resonance scattering experiment is already determined by the pole term (|2.8|) . 
and is therefore given by 

w r . T . = 2xr£( n ; i ) ( - i r^t©i4f )( "-"< + ^i (3.9) 

n=0 ^ ' ' k=0 ^ ' 

n=0 ^ ' 



up to a normalization factor which will have to be determined by normalizing the overall 
probability to 1. Here we define the operator 

^. = ^Eft)i4) (fc)(n - fe) ( + ^|. (3.10) 

A;=0 ^ ' 

We hypothesize that the microphysical state from higher order poles connected with the 
decay experiment has the same structure as the microphysical state (|3.9| ), which is certainly 
in agreement with the first order case ( p . 8| a) in comparison with (3.8b), 



w = L r +1 ) <-o-£e (i) m i) "-"t^i on) 

n=0 ^ ' ' k=0 ^ ' 



For r = 1, this reduces to (|3.8| a). It should be mentioned that mathematically there is an 
important difference between ( |3.9|) and ( |3.11 ) because the {i[)~\z 1 )^ ^ n ~ k \ + z\<p + ) are analytic 
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functions for z in the lower half-plane, whereas the (?/>f \z~f k ^ ( n ~ k \~ z\ip2) are n °t- 
Whether the microphysical state of the (hypothetical) quasistationary microphysical system 
is always represented by the mathematical object ( |3.11| ) or whether also each individual 



W(n) = V ^X{i) l^) (fe) (n ~ k) (~ Z «\ : - = 0, 1, . . . , r - 1 , (3.12) 

has a separate physical meaning, cannot be said at this point. 

This means that the conjectural physical state associated with the r-th order pole is 
a mixed state W, all of whose components W^ n \ except for the zeroth component W^°\ 
cannot be reduced further into "pure" states given by dyadic products like \z^f k ' ^ k \~zji\. 
This is quite consistent with our earlier remark that the label k is not a quantum number 
connected with an observable (like the suppressed labels 6 2 ,...,6 n ). Therefore a "pure 
state" with a definite value of k, like \z R ^f k " > ^\~Zr\, k > 1, does not make sense physically. 
A physical interpretation could only be given to the whole r-dimensional space M. ZR , ( |2.13| ). 
The individual W^ n \ n = 0, 1, 2,...,r — 1, act in the subspaces M.zr C -M Zr which 
are spanned by Gamow vectors of order 0, 1, . . . ,n (Jordan vectors of degree n + 1, i.e. 
(H x — Zfj) n+1 A4 z n J = 0). Here the question is whether there could be a physical meaning to 
each separately, or whether only the particular mixture W given by ( |3.11| ) can occur 
physically. 

Though the quantities |^)^ ^ k \~zji\ will have no physical meaning, even if higher order 
poles exist, they have been considered in the literature (Antoniou and Gadella, 1995; Golod- 
berger and Watson; Newton, 1982; Goldhaber, 1968), and their time evolution is calculated 
in a straightforward way: 



(e mt y\z R f^ k \-z R \e mt = (3.13) 

■"E E (f) (*) Htym m \z R f- l)(k - m \-z R \ . 

;_n ™— n V / V / 



1=0 m=0 



This state operator shows the additional polynomial time dependence that has always been 
considered an obstacle to the use of higher order poles for quasistationary states. A polyno- 
mial time dependence should have shown up in many experiments. 



We now derive the time evolution of the microphysical state operators ( 3.11| ) from higher 



order poles of the S-matrix using the time evolution obtained for the Gamow- Jordan vector. 
It will turn out that the operator ( |3.12| ) and therewith ( |3.11| ) have a purely exponential time 
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evolution. Inserting ( |2.22| ) and ( 2.2'd ) into 

WW(t) = (e iHt ) x W^e iHt = (fj {e tH T Kf k) (n-fc) t*Jil <? H \ t>0, (3.14) 

fc=0 ^ ' 

we get 

^(t)=e-v^^Wj (*) ( n ~ k ) {-rtr i ^- k --\ ZR r^\- ZR \ . 

k=0 1=0 m=0 V / V / V / 

After reordering the summations and the terms in the binomial coefficients one can separate 
out the dyades, 

-pn n / \ n ~ m / \ n - m / A 

w(n) (*) = ^ E u) E n j lz « fl) (m)( ~^ S [Tl) H*) fc -^r fe ~ m • 

m =o ^ ' 1=0 ^ ' k=l ^ ' 

Since the indices labeling the Gamow- Jordan vectors do not depend upon k, the sum over k 
may be performed using the binomial formula 

E/n-m-A ,. f }k-i ( if \n-k-m _ t if _ i + \n-m-i _ f 1 for Z = n - m 1 
^ Jfe-Z y 1 J 1 J ~ l J ~\0 for l^n-m J ~ 

and one gets 

■pn n / \ 

= e - rt — E l^) ( " _m) M t^l = e" r * W^ (n) (0) ; t > . (3.15) 



m=0 x ' 

This means that the non-reducible (i.e. "mixed") state operator of (|3.12| ) has a simple 
exponential semigroup time evolution, and also the operator W of ( |3 . 1 1| ) being a linear 
combination of the W^ n \ 

W(t) = (e iHt ) x We lHt = e- Tt W , t > . (3.16) 

It turns out that the operator fl3.12| ) is the only operator in Mj^ formed by the dyadic 
products \zr) ^ m '^\~zji\ with m, I = 0, 1, ••• ,n, which has a purely exponential time 
evolution, thus being distinguished from all other operators in A4%£. Thus we have seen 
that the state operator which we conjecture from the r-th order pole term describes a non- 
reducible "mixed" microphysical decaying state which obeys an exact exponential decay 
law. 

In analogy to (|3.6j ) one can also calculate the time evolution of the operators (|3.9| ) 
and ( |3.10| ), but their time evolution will always have an additional polynomial time de- 
pendence besides the exponential. 
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4 General Form of the Exponentially Decaying Gamow 
State 

In this section we want to discuss the converse of the above reasoning, where we conjectured 
the density operator for higher order decaying Gamow states and derived a purely exponential 
time evolution. We ask the question: If we require exponentially decaying time evolution 
for a Gamow state operator formed by dyadic products of vectors in A4 ZR , what is the most 
general form of such an operator? 

We denote the most general form of W of finite dimension j e IV by 

3 3 
fc=0 h=0 

with arbitrary coefficients A^k- We want to change the order of summation from the states 
of the form \z^f k ' ^ h '{~ Zr\ to the states of the form \z^f k ' ^ n ~ k \~ZR\. Changing the label 
h = n — k we write this sum as 

3 3+k 
k=0 n=k 

Switching the order of k and n divides Q4.2| ) into two sums, 

j j+k j n 2j j 

= E E • («) 

k=0 n=k n=0 k=0 n=j+l k=n—j 

We define = + W { f )} such that 

K) = tt^|%f )M r^|, (4.4) 

n=0 k=0 

K) = EE A n-kAz- R ) [k) {n -\z R \ . (4.5) 

n=j+l k=n—j 

In the following, we calculate the coefficients of and give an argument why this suffices 
to conlude that the coefficients of Wf~{ all turn out to be zero. 
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The time dependence of is given, using ( |2.22| ) and ( |2.23| ), by 

3 n 
n=0 k=0 

j n k n—k / . \ / . \ 

= e " r< EEEE A ^ (i) r~ k ) t-*r i ^r k ~ m \^ R ) {1) (w) t^i 



n=0 fc=0 i=0 m=0 

Changing the order of the summations, 

jnk n—k j n n n—k j j n n—k j j n—l n—m j j—l j n—m 

eeee=eeee=eeee=eeee=ee e e. 

n=0 fc=0 (=0 m=0 n=0 (=0 fc=Z m=0 «=0 n=i fc=Z m=0 1=0 n=l m=0 fc=Z (=0 m=0 n=l+m k=l 

allows the dyadic products, which are linearly independent operators, to be factored out of 
the sums over terms in which they appear as common factors: 

j j-l j n—m 



w &® = e " rt EE E E A ^ (?) ( n ~ k ) w n - k - m \z- R p {m \-z R 

1=0 m=0 n=l+m k=l V / V / 

3 j- 1 3 n-m /l \ / 7 \ 

= e - rt EEi^) (z)(w) t^i E { u (^r'"'" 

Z=0 m=0 n=«+m fc=Z V / V / 

J J — i J n—m /t \ / 

= e - rt EEi%> (z)M r^i E (^r m "E^(")( 

Z=0 m=0 n=l+m k=l \ / \ 



A; \ / 72 — A; 

m 



The operator (t) will decay according to the pure exponential e~ rt if and only if all 
terms involving additional powers of t cancel. All terms involving additional powers of t will 
cancel if and only if the coefficients A n _^^ satisfy the conditions 

I e {Or- • 

me{0r-- J-l-l}, (4.6) 



fc=; \ / \ / [ fie{ra + ! + l,---,j}, 

The simplest of these conditions are those for which n = m + I + 1, i.e., those for which 
m — n — I — 1, because they are the only conditions that involve sums over only two values 
of k: 

1+1 



A "- l > 1 \n-l-l) ~ A "- l - 1 ' l + 1 { I ) tor { I e {0, ■ ■ ■ ,n - 1} 
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or, equivalent ly, 

or, equivalent ly, 
A 

- rL n—K,K 



n — k + 1 
n — k 



A 



n—k,k 



k 
k-l 



for 



7i-AH-l)!(Jfe-l)! 



(n-k)\k\ 



A, 



n—k+l,k—l 



for 



n e {1, ■ ■ ■ ,j}, 
k e {!,-•• ,n} 



n g {1,-- • , j}, 
fc G {!,■■• ,n}. 



These conditions relate pairs of coefficients A n _k,k having the same values of n and successive 
values of k. For fixed n G {1, 2, • • • , j}, they may be used recursively to show that A n _k,k 
must equal A Hj o multiplied by the binomial coefficient (^) = k ^i k y/ - 



A 



n—k,k 



(n-k+l)\(k-l)\ 



(n-k)\k\ 



(n-fc+2)!(fc-2)! 
}n-k+l)\(k-l)\_ 



nW 



{n-k)\k\ 



A 



n,0 



n 



A 



1 1,0 



for 



(n-l)!l! 
_(n-2)!2!_ 

n g {l,- • • , j} 
k e {!,■•• ,n} 



n\0\ 



_(n-l)!l! 



A n ,o 
(4.7) 



Substituting this result into the full set of conditions Q4.6Q , using the identity 



n\ fk\ fn — k 



k \l 



rn 



n\ fn — m\ fn—m — l 



rn 



I 



k-l 



and then using the binomial formula gives 



° - ^"g(:)C)( 



n\ fk\ fn — k 

m 



-1) 



k-l 



A n ,o 

A n ,o 
A n ,o 
A n ,o 



n \ n—m 



rn 



n \ / n—m 



rn 



n \ n—m 



E 

k=l 
n—m— I 

E 

k-l=0 



n—m — I 
k-l 

n — m — I 
k-l 



-1) 



k-l 



I n—m—k 



1) 



k-l 



rn 



n \ I n—m 



(1-1)" 
n for 



/e{o,---,j-i}, 

m G {0,-- - ,j - 1 -/}, 
n G {m + / + !,••• ,j}, 
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which shows that the remaining conditions are automatically satisfied by (|4.7|) without plac- 
ing any further conditions on the coefficients A n $. The coefficients for n £ {1, • • • , j}, 
and also the coefficient Aq^, remain completely arbitrary. 

We conclude that a linear combination of dyadic products \zr)® ^ m \~ZR\ decays according 
to the pure exponential e~ Tt if and only if it is of the form 

n=0 k=0 ^ ' 

with arbitrary coefficients A n $- 

Now, coming back to the argument at the beginning of this section: We have shown that 
for arbitrary j, the operator WB\ depends only on the choice of the coefficients v4 n>0 . Since 

j was chosen arbitrarily, we can also take an operator W of dimension 2j, such that it 
contains at least all the terms belonging to the operator of ( f4.1|) and some additional 
terms which we set zero by choice of the coefficients A^ = for h,k > j, 

W§) = W^ j} with A hk = for k > j or h > j (4.9) 

2j 2j-k 

= Yl A hk\ z R) {k) {h X~ z R\ ; with A hk = for k > j or h > j . 

k=0 h=0 

If for n > j the A n Q = then we know that, according to ( |4.7| ), all the other terms A n _k t k 
are zero. Since they make up the the operator W^l of (|4.5|) , we conclude that the coefficients 

of are zero, as demanded above. 

Comparing ( |4.8|) with arbitrary coefficients A n fi to the Gamow state operator ( |3.11| ) 
suggested by the pole term, one sees that the structure is the same with j = r — 1 and 

4M. = a * r („ + i)H)*£- MO) 

5 Conclusion 

Gamow vectors that can describe resonances and decaying states from first order poles of the 
S-matrix have been known for two decades. In this paper, we discussed their generalization 
to Gamow vectors describing resonances from higher order S-matrix poles. This led to a set 
of r higher order Gamow vectors associated to a pole of multiplicity r, which are Jordan 
vectors to a self-adjoint Hamiltonian with complex eigenvalue Er — iT/2. They are basis 
elements of a generalized eigenvector expansion which suggests the form of a microphysical 
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state associated with this higher order resonance pole. This microphysical state is a mixture 
of non-reducible components, and in spite of the fact that the higher order Gamow vectors 
have an additional polynomial time dependence, this microphysical state obeys a purely 
exponential decay law. We showed that this state operator has the same structure as the 
most general from constructed from decaying higher order Gamow vectors that leads to an 
exponential decay law. 

Katznelson, 1980; Bhamathi and Sudarshan, 1995; Brandas and Dreismann, 1987; Antoniou 
and Tasaki, 1993). It has been shown that Jordan blocks arise naturally from higher order 
S-matrix poles and represent a self-adjoint Hamiltonian by a complex matrix in a finite 
dimensional subspace contained in the rigged Hilbert space. Although higher order S-matrix 
poles are not excluded theoretically, there has been so far very little experimental evidence 
for their existence, because they were always believed to have polynomial time dependence. 
Our results suggest that the empirical objection to the existence of higher order poles of the 
S-matrix does not rule out the possibility of exponentially decaying states constructed from 
higher order Gamow vectors. 
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